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Abstract 

We show that the von Neumann's algorithm of reduction (i.e. the algorithm of calculating the 
density matrix of the observable subsystem from the density matrix of the closed quantum system) 
corresponds to the special approximation at which the unobservable subsystem is supposed to be 
in the steady state of minimum information (infinite temperature) . We formulate the generalized 
algorithm of reduction that includes as limiting cases the von Neumann's reduction and the self- 
congruent correlated reduction most corresponding to the quantum nondemolition measurement. 
We demonstrate the correlation in dynamics of subsystems with exactly soluble models of quantum 
optics: 1) about the dynamics of a pair of interacting two-level atoms, and 2) about the dynamics 
of a two-level atom interacting with a single-mode resonant field. 
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I. INTRODUCTION 



The old question of the quantum theory on the reduction of a quantum system and on 
the decoherence of an observable subsystem (transition from a pure state to a mixed state) 
due to its interaction with an unobservable subsystem and the measurement is formulated 
anew in quantum optics [1, 2]. This process is initiated by problems of quantum comput- 
ing [3], by hope to clarify the quantum nature of the irreversibility [2, 4] and, mainly, by 
the new experimental potentials to check consequences of axioms of quantum mechanics 
accepted previously. The qualitatively new experiments allow one to observe dynamics of 
subsystems of simple quantum systems (e.g., the dynamics of a pair of correlated photons 
emitted in a nonlinear process, or the dynamics of a single two-level atom in the course of 
its interaction with the resonant field of extremely small intensity [1, 2, 3, 5]). The math- 
ematical operation of reduction defines the rule of calculating the density operator of an 
observable subsystem from the known density operator of a closed quantum system. By the 
von Neumann's postulate, this rule consists in calculating the partial trace of the density 
matrix of a closed quantum system over the unobservable subsystem [6, 7, 8, 9]. Though 
this nonunitary operation has been universally accepted to account for the influence of the 
classical measuring device on the result of measurement, to realize the transition from the 
closed system to the open one, and to introduce the irreversibility into the dynamics of 
an observable subsystem [2, 4, 6, 7, 9], the algorithm of accounting for the interaction of 
the closed quantum system with a measuring device remains in so doing reduction implicit. 
Moreover, when calculating the mean values of observables of subsystems, the realization of 
reduction is imphcit ("disguised") by itself. 

The initial point of the our treatment lies in the identity of the mathematical formulation 
of the reduction (other than the interpretation!) with the formulation of the quantum 
relaxation theory problem [8, 9]. In both cases the closed composite quantum system is 
considered. Consequently the unitary time evolution of the system is described by the 
density operator p{t), which may be approximately calculated (basically in general case). 
In both cases the dynamics, i.e. the density operator, of one subsystem Pait) of the system 
should be calculated. This extraction of the one subsystem implies implicitly the separation 
of the system into two parts: the observable subsystem a and the unobservable subsystem 

In general case the subsystems of the arbitrary closed system are equal in rights. Let 
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us suppose the dynamics of the subsystem (5 has described by some density operator P0{t). 
The exact representation of the density operator of the closed system through the density 
operators of subsystems is known for the noninteracting subsystems only: p{t) = Pa(t)®Pf3(t) 
[8, 9]. The considered problem consist in calculating of the density operator of one subsystem 
in the general case of interacting or interacted previously (entangled) subsystems. 

In the quantum relaxation theory the " unobservable" subsystem is believed to be the 
thermostat at a temperature of T from physical reasons. Consequently its density operator 
supposed to be known and stationary: ppiT). The desired density operator Pa{t) next is 
calculated to some approximation from the approximate equation p{t) ~ paif) ® PpiT), 
where two operators /3(t) and PpiT) are known. In the reduction problem the state of 
the unobservable subsystem is believed unknown. The algorithm was instead postulated 
to calculate Pa{t) from p{t). The question is arises then: is this postulation equivalent to 
some approach on the state of unobservable subsystem? As a matter of fact this question is 
identical with the mathematical problem of the quantum relaxation theory: how the density 
operator of subsystem Pi3{t) should be approximately calculated from the known density 
operators p(t),pQ,(t)? 

We apply the method similar to used by von Neumann for demonstrating the correlation 
(the one-to-one correspondence) between the physically observables of the different inter- 
acting subsystems of the closed quantum system [6] . In this way the approximate equation 
system demonstrating the mutual relations between the density operators of subsystems and 
the density operator of the closed system was deduced. This equation system reveals that 
given in some way the density operator of one subsystem, the density operator of the other 
subsystem can be calculated approximately. We prove that the von Neumann's algorithm 
of reduction is identical to the approximation when the subsystem /3 is constantly in the 
state of minimum information ppmin (identical with the state of infinite temperature [6, 8]). 
It should emphasize this approximation is solely employing when the reduction earring out, 
other than when the p{t) calculating! 

By this means the von Neumann's reduction is mathematically equivalent to the replace- 
ment p{t) Pa{t) ® Pf3min with cxactly defined operators Pa{t)^ Pfimin- This replacement 
corresponds on other hand to the changing of the considered quantum system and to the 
traditional interpretation of reduction as the transition from the closed to an open system 
or the accounting for the measurement. The state of minimum information (of infinite tem- 
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perature) corresponds to the complete absence of the quantum coherence in the subsystem 
13. This imphcit approach leads to some decoherence of the observable subsystem although 
all the calculations are done in framework of the closed quantum system. 

We get the conclusion that any one of the possible algorithms of reduction is mathemat- 
ically equivalent to the replacement f){t) — >■ pQ,{t) ® Pjiif)-, to the approximate represen- 
tation of the closed quantum system as two quasy-independent subsystems with the definite 
dynamics. This raises the question of whether the accepted algorithm of reduction (i.e. the 
corresponding replacement of the density operator) corresponds to the initial assumption of 
the closed quantum system. Or mathematically: how much would this replacement be in 
error with the believed to be true p{f)l 

As it is evident from ours consideration, the generally accepted von Neumann's algorithm 
of reduction is self-contradictory: on the one hand the closed quantum system is assumed, 
and on the other the infinite temperature of the subsystem is taken implisitly. It is possible 
to perfect the algorithm of reduction: for example, to accept some more reasonable state of 
the unobservable subsystem (the thermostat, as in the quantum relaxation theory, or one of 
eigenstates, as in the theory of quantum nondemolition measurements of photons [12, 13], 
or any intermediate variant). But in any case this step corresponds to the changing of the 
considered system and to the coming to the new system with the new dynamics, what is 
customary interpreted as the coming to the open system. 

If the both subsystems of the considered system are the essentially quantum systems and 
the approximation of the closed system is believed true, it seems reasonable to use the prin- 
ciple of the interchangeability of subsystems [6] to make the reduction. We use the mutual 
relations between the density operators p(t), pa{t), pfs(t) to calculate the correlated density 
operators pacit)i Pfscif) simultaneously by the self-congruent procedure of the successive 
approximations. In this way the reduction making the minimum (the vanishing in the limit) 
changing of p{t) is obtainable. In other words, the correlated reduction procedure is most 
consistent with the initial approximation of the closed quantum system or to the quantum 
nondemolition measurement. 

The structure of the paper is as follows: in Sec. II we demonstrate that the von Neu- 
mann's reduction corresponds to the approximation of the steady state of minimum infor- 
mation of the unobservable subsystem and formulate the generalized reduction algorithm 
as an approximation of a given state of the unobservable subsystem of a closed quantum 
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system. Sec. Ill contains the algorithm of correlated reduction that defines the procedure 
of calculating self-congruent density operators of subsystems of a closed quantum system 
in the limit of the quantum nondemolition measurement. In Sec. IV we use the correlated 
reduction technique to consider three known problems of quantum optics: 1) the separable 
representation of the entangled EPR(Einshtein-Podolsky-Rosen) state, 2) the dynamics of a 
pair of interacting two-level atoms, 3) the dynamics of a two-level atom interacting with the 
resonant single- mode electromagnetic field (the model of cavity quantum electrodynamics). 
Designations and definitions [6, 7, 8, 9]: 

1. H — Ha -\- Hp -\- Hce/3 is a hamiltonian of a closed quantum system. Ha and Hp 

are hamiltonians of uncoupled subsystems a and i3, respectively, Hap is the interaction 
hamiltonian of subsystems; A and B are operators of observables of subsystems a and (3, 
respectively. 

2. p{t) — U{t)p{0)U^t) is the solution of the Schrodinger equation determining the 
dynamics of the system H, U{t) — e:K.p[—{it/h)H] is the operator of evolution. 

3. Operators Ha, A and Hp^B are determined in spaces of subsystems a and /3, respec- 
tively; operators H, Ha/s, p{t),U{t) are determined in the complete space that is the direct 
(exterior ^ (g)) product of spaces a, (3; A' = A ^ ^p,B' = la <^ B are extended operators 
of subsystems in the complete space; la{i3) are unity operators in spaces of subsystems, 
la (E) 1/3 = 1 is the unity operator of the complete space; a^, bi'f are matrix elements of the 
operators A, B; o^u'^jfy are matrix elements of operators (O) of the closed system, a pair 
of subscripts {ii'} designates basis state of the closed system. 

4. The von Neumann's algorithm of calculating the reduced density operator of the 
observable subsystem is determined by: 

PaNit) = Spf3p{t), or Pl3N{t) = Spap{t), (1.1) 

where Spa[i3) means the partial trace over the subsystem q;(/3), the subscript N hereinafter 
designates quantities calculated in accordance with the von Neumann's algorithm of reduc- 
tion. 

II. THE THEOREM ON THE CORRELATION OF SUBSYSTEMS 
Theorem 
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The reduction (1.1) of the closed quantum system H for the observable subsystem a cor- 
responds to the approximation at which: i) the closed system is represented as two uncoupled 
correlated subsystems, and ii) the unobservable subsystem (3 is supposed to be in the steady 
state of minimum information. 

Proof. 

i) . Let us suppose that the dynamics of the closed system H is approximately represented 
as the dynamics of two uncoupled subsystems. It means [8, 9] that there are such density 
operators Pac{t) and ppc{t) {SpaPac{t) — Sppppc{t) — 1) of subsystems for which 

Pac{t)®hc{t)^p{t)- (2.1) 

We transform (2.1) in much the same way as von Neumann did to prove the mutual corre- 
lation of mean values (^(i)), {B{t)) ( see [6, Ch.6.2]). If we multiply (2.1) on both sides by 
Ppci^) — ^cc® Pf3c{t) or Pac{t) — PaC <8) 1/3 ^ud take partial traces over the subsystem /3 or 
a we obtain the system of two connected equations [10, 11]: 

SPapp{t)ppcV') 
^ /,N _ Sp^p{t)p^cit) o^ 

Ppc{t) ^ ^ ^, . 2.3 

Spal3p[t)p^C\t) 

The right-hand sides of (2. 2), (2. 3) represent the normalized quantum averaging of the density 
operator of the closed system over one of the subsystems. Each equation determines the 
density operator of one subsystem through the density operator of the closed system and 
density operator of another subsystem. Prom (2.1) and (2. 2), (2. 3), it is evident that density 
operators Pac{t), P/scit) are correlated, that is why we may call i) the approximation of 
correlated subsystems (is marked by a subscript C). 

ii) . Let us suppose that the unobservable subsystem j3 is in the steady state of minimum 
information, i.e. its density operator equals [6, 8]: 

hc{t) = hmin = (l/iV/3) (2-4) 

where Np is the dimension (number of eigenstates) of the subsystem /3. The theorem then 
follows if we insert (2.4) into the right-hand side of (2.2): 

Pa{t) ^ ^ ;V ^ = Sp^m = PaN{t). (2.5) 



Sp^f3p{t){i/Nf^) (ia^ip) 



Consequence 1 

By (2.1)-(2.4), the von Neumann's algorithm of reduction (1.1) is equivalent to the re- 
placement of the density operator p{t) with the reduced density operator pRN{t) according 
to the rule: 

p{t) ^ PRN{t) = PaNit) ® PfSmin = [Sppp{t)] ® (V^/j) (2-6) 

The relation (2.6) agrees in form with " the basic equation of irreversibihty " or "the first 
assumption" of the quantum relaxation theory ( see, e.g. [9, Ch.7]): 



exp 




{l/kBT)Hf, 




Spi3 exp 


'-{l/kBT)Hp 



that redefines p{t) in a case, when the unobservable subsystem is taken to be the thermo- 
stat at a temperature of T {ks is the Boltzmann constant). It follows from the comparison 
of (2.6) with (2.7) that the approximation ii) within the quantum relaxation theory corre- 
sponds to the representation of the unobservable subsystem by the thermostat at the infinite 
temperature, no matter of actual physical properties of the unobservable subsystem. 
Consequence 2 

i), and ii) represent qualitatively distinct approximations. The first approximation (i) 
represents a closed quantum system as two of such quasi-uncoupled correlated subsystems 
with mean values of all observables of the closed system being approximately correct. The 
second approximation (ii) describes the action of the measurement on the state of the ob- 
servable subsystem and implicitly carries out the transition from a closed system to an open 
one. By Consequence 1, the model of measurement (1.1)=(2.6) used in the von Neumann's 
reduction may be considered as the identification of the classical measuring device with a 
thermostat of the infinite temperature and the measurement with the ideal thermal contact 
of the unobservable subsystem with this device. The limited applicability of such model and 
hence the von Neumann's reduction seems physically evident. 

In the form similar to ii), the model of measurement was included into the algorithm of 
reduction in the theory of quantum nondemolition measurements of quantum optics [12, 13]. 
The "probe" (unobservable) subsystem (3 was postulated to be determined with certainty 
in its eigenstate at the measurement point of time t, i.e. the probe subsystem density 
operator is to be presented by a projector ppj'it) — Ppj' = The equation of such 
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"projective" reduction (is marked by the subscript P) was determined by [12, 13]: 

being the special case of the equation (3a). 

Uniting (2.6) = (1.1) and (2.8), it is possible to formulate the generalized algorithm of 
reduction as the procedure of calculating the density operator of an observable subsystem 
in the approximation of a given state of the unobservable subsystem in the form of two 
consecutive operations: 

1. A model of interaction of a closed quantum system with a measuring device 
is postulated from physical reasons and is expressed mathematically in the 
form of the approximation of a given state of an unobservable subsystem. 

2. The density operator of the observable subsystem is calculated from equa- 
tion (2.2) of correlated quasi-uncoupled subsystems of a closed quantum 
system. 

In (2.6) and (2.8) two special approximations of the given state of the subsystem (3 
are realized. In (2.6) the subsystem (3 is permanently in each of eigenstates with equal 
probability and in (2.8) it is with certainty in one of them at the time t. 

Consequence 3 

It is known [8, 9] that calculating the mean value of the observable of a closed system 
described by the operator A', the von Neumann's reduction of a closed system is carried out 
implicitly: 

(i'(t)) = Sp^pp{t)A = Spp^N{t)A = (i(t))iv. (2.9) 

In doing so not the state of the subsystem /3 is postulated but the operator of the closed 
system that corresponds to the observable A of the independent subsystem a {A' -^-^ A) (see 
[6, Ch. 6.2]). Hence, by Theorem and Consequence 2, the choice of the required operator 
as A' is equivalent to realization of both approximations i) and ii). 

Such equivalence becomes evident by considering the correlator of observables A and B 
of subsystems (^(^A^ (i)^ = Spa/3p{t) (^A(S> B^- If we use the obvious identity A(S> B — 
(^A^il3^ (ia^B^, the commutativity A',B' — and the interchangeability of partial 
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traces SpajS — Sp^a, the correlator can be presented in two forms symmetric relative to 
subsystems: 



(2.10) 



where 



A{t))^ = SpaPaB{t)A, {B{t))^ = SppppA{t)B, (2.11) 

SpaPp{t)PpB Spapp{t)p'^A 

Po^A^A/Spo^A, ppB^B/SpfsB, (2.13) 



It is supposed here for simplicity that (^A{t)j^ ^ 0, (^B{t)J^ ^ 0, and ^4, B are the nonneg- 
ative operators. 

The expressions (2.10)-(2.13) are the exact equations for arbitrary Hermitian ^4,5. 



The meaning of values {^Aityj ^ and {Bityj ^ is clear from the comparison of (2.12) with 
(2. 2), (2. 3). For example, the operator paB(t) is the reduced density operator of the sub- 
system a in the approximation of the given state p^B of the subsystem (3 in accordance 
with Consequence 2. Accordingly, (^A{t)'^^ is the mean value of the observable A of the 
subsystem a if the subsystem (3 is in the state p^^. Thus, (2.10) shows that the correlator 
(^(A0 B^ (i)^ of the closed system H is equal to the product of mean values of observables 
of subsystems calculated for the system H reduced in the different ways. Namely: 1) the 
reduced state of a subsystem a is defined at the given = pas, but the reduced state of 
the subsystem f3 is defined simultaneously at the given p^ = Pamm, or, leading to the same 
result, 2) the reduced state of the subsystem a is defined at the given p^ = Pjjmin but that 
of the subsystem (3 at pa = PaA- 



In the specific case S = l^j we have p^B = P/Smin: (-^(0/^^ = \^(*^/a ^ \^^^\ b 
{A{t))^ and (2.10) = (2.9). 

For this axiomatic choice of the operator of an observable of the closed system A' <-> 
A, commonly accepted in the quantum mechanics, in the limits of uncoupled subsystems 
{Hafi = 0) the expression of the observable {A'{t)) = {A(t)) is exact. But at Hajs 7^ the 
mathematically exact equality (2.9) is the physical approximation as it contains implicitly 
approximations i) and ii). One can see from (2.6), (2.8) - (2.13) and Consequence 2 that 
depending on a choice of a model of measurements (or of an algorithm of reduction) the 
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calculated mean value of the observable described by the same operator A and the same 
density operator of the closed system p{t) is different: {A{t))N 7^ {A{t))pj' ^ {A{ty)B- 
The physical reason of this difference is obvious — a choice of a model of measurement 
(of a algorithm of reduction) corresponds to a choice of the back effect of the measuring 
device on the state of the observable subsystem and, hence, on the result of measurement. 
Or, alternatively, it is possible to say following (2.10) that each choice of the model of 
measurement corresponds to the calculation of different observables of the closed quantum 
system. 



III. ALGORITHM OF THE CORRELATED REDUCTION 



If the approximation of the closed system is true, it seems reasonable to seek the algorithm 
of reduction such that minimize the changing of p(i). 
Let us consider expressions: 



m)c = sppacm, {m)c = sphcim (3.1) 

((i S) {t))^ = [Spp^c{t)A\ [Spppc{t)B] . (3.2) 

The equalities (3.1) define mean values of observables A and B of independent subsystems 
with the correlated density operators Pac{t) and ppc{t), respectively. But these correlated 
density operators of subsystems can be determined only approximately for the lack of the 
exact expression of p{t) in the form of (2.1) at Hap 0- Let us calculate Pacit), P/scit) from 
the system of equations (2. 2), (2. 3) by the successive approximation method [10, 11]: 

Sppp{t)pf\t) 

,(n+l),,, _ Sp^p{t)p^:\t) 
Spat3p{t)pa\t) 

where, for an example, the von Neumann's reduced density operators PaNif)-, Pi3N{t) may 
be used as a zeroth-order (n=0) iteration. For the convergent sequences (3. 3), (3. 4), the 
calculated self-congruent density operators of subsystems most closely correspond to (2.1). 
It means that the mathematical operation of reduction (3. 3), (3.4) does not change approxi- 
mately the state of the closed quantum system and hence may be considered to correspond 
to the quantum nondemolition measurement. 



P<^c{t) = jim + \t) = jim ITTTHTTT' (3-3) 



- r -(n+l /,x ^PaPyi^ Pa V^ roA\ 
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A matter of principle is if (2.9) or (3.1) describes observables of subsystems of a closed 
system more properly. 

The expression (2.9) is the mathematical equality but the physical approximation by 
Consequences 1-3. Each of (3.1) is the physical equality but mathematical approximation by 
(2.1) and (3.3), (3.4). The accuracy of (3.1) or (3. 3), (3. 4) can be estimated by a comparison 
of matrix elements at the left-hand and right-hand sides of (2.1) or by a comparison of mean 
values at the left-hand and right-hand sides of (3.2). It should be particularly emphasized 
that the correlated reduction (3. 3), (3.4) seems to be more appropriate to the equivalency of 
subsystems of the closed system than the reduction (1.1) including the physically impossible 
property (T = oo) of the unobservable subsystem. Therefore we expect that the correlated 
self-congruent density operators (3.3), (3.4) would more properly describe the dynamics of 
subsystems in experimental cases when the approximation of the closed system holds. 

Let us consider the previous procedure of a correlated reduction and its physical 
consequences on model problems of quantum optics. 

IV. DYNAMICS OF A PAIR OF CORRELATED TWO-LEVEL SUBSYSTEMS 

Two coupled two-level subsystems constitute the compound closed quantum system now 
in use to describe many basic problems of quantum optics and quantum informatics [1, 2, 3]. 
Let us look after the dynamics of such simplest correlated subsystems. To visualize the 
distinction in the dynamics of subsystems described by the traditional reduction (1.1), (2. 6), 
and by the correlated reduction (2. 2), (2. 3), (3. 3), (3. 4) we present at first the consideration 
in general. 

Let us designate cigcnstates of uncoupled subsystems a and /5 by |2)q,, |1)q, and |2)^, |1)^, 
respectively. Operators of the closed system are represented by matrices of the fourth order, 
each element of which is designated by two pairs of subscripts (see, e.g., [9, Appendix A]). 
The first subscript of each pair is attributed to a subsystem a and second subsctipt to a 
subsystem /3. The general expression of the density matrix of the considered closed system 
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is then: 



^ P22,22 P22,21 P22,12 P22,ll ^ 



(4.1) 



P21,22 P21,21 P21,12 P21,ll 
Pl2,22 Pl2,21 Pl2,12 Pl2,ll 

y Pii,22 pii,2i Pii,i2 pii,ii y 

For an observable subsystem a the reduced von Neumann's density matrix in accordance 
with (1.1) is defined by expression: 



[paN) 



P22,22 + P21,21 P22,12 + P21,ll 
Pl2,22 + Pll,21 Pl2,12 + Pll,ll 

Accordingly, for an observable subsystem /3: 



(4.2) 



^ P22,22 + Pl2,12 P22,21 + Pl2,ll ^ 



(4.3) 



P21,22 + Pll,12 P21,21 + Pll,ll y 

In accordance with (2.6), such definition of the operation of reduction is equivalent to the 
replacement of (4.1) by expressions: 



(P) 



{PBxxn) — {Po.n) ® {PfSmin) 
( 



P22,22 + P21,21 P22,12 + P21,ll 

P22,22 + P21,21 P22,12 + P21,ll 

Pl2,22 + Pll,21 Pl2,12 + Pll,ll 

Pl2,22 + Pll,21 Pl2,12 + Pll,ll j 



(4.4) 



or 



(P) 



{PR^n) = {Pamin) ® (PPn) 
I 



\ 



(4.5) 



P22,22 + Pl2,12 P22,21 + Pl2,ll 

P21,22 + Pll,12 P21,21 + Pll,ll 

P22,22 + Pl2,12 P22,21 + Pl2,ll 

P21,22 + Pll,12 P21,21 + Pll,ll I 

For comparison, we write the system of equations (2. 2), (2. 3) determining correlated reduced 
density matrices of subsystems of the closed system (4.1): 



V 



(Pac) 



' a22 0(21 ^ 

^ (Xi2 ail J 



^ (P22, 22/^22 + P21,2l/5ll+ (p22, 12/322 + P21,ll/5ll+ ^ 

+P22,2l/?12 + P21,22/?2l) +P22,ll/?12 + P21,12/32l) 

(Pl2,22/522 + Pll,2l/5ll+ (Pl2,12/322 + Pll,ll/3ll + 

+Pl2,2l/5l2 + Pll,22/32l) +Pl2,ll/3l2 + Pll, 12/^21) y 
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X (^22,22/522 + P21,2l/3ll + P22,2l/3l2 + ^21,22/321 
+Pl2, 12/322 + + Pl2,llA2 + Pll,12/32l) ^ , 



(4.6) 



1322 1321 
/5l2 Al 



(P22,22Q;22 + P12,12Q;ii+ (P22,21Q;22 + Pl2,llQ;il + 
+P22,12ttl2 + Pl2,22a2l) +P22,liai2 + Pl2,2ia2l) 
(P21,22«22 + Pll,12ail+ (P2L22a22 + Pll,liail + 
+P21,12ai2 + Pll,22Q!2l) +P21,llQ!l2 + Pll,2lQ!2l) y 
X (P22,22Q;22 + P12,12Q;ii + P22,120ll2 + Pl2,22Q;21 

+ P21,21Q:22 + Pi1,11Q:11 + P21,llQ:i2 + P\l,2l0.2\)~^ ■ 



(4.7) 



It is seen that in the general case the mutual correlation of subsystems (4.6) and (4.7) as 
distinct from (4.2) and (4.3) implies the following. 

1) The possibility of the simultaneous existence of quantum coherence, i.e. nonzero nondi- 
agonal elements of reduced density matrices of both subsystems (q;i2 7^ 0, /3i2 7^ 0). 

2) The inequality of probabilities to find the subsystems in their eigenstates (an ^ 0:22, 7^ 

(322). 

Let us use (4.1)-(4.7) in considering two known physical problems. 



A. Representation of an entangled EPR (Einstein-Podolsky-Rosen) pair of pho- 
tons as the superposition of separated correlated photons 



The model of the entangled EPR state is of considerable current use in the description 
of experiments on the correlation of photons [1, 2, 3]. The density matrix of such closed 
quantum system of two noninteracting but correlated photons in designations (4.1) is of the 
form [1, 2, 3]: 



'^0 0^ 



1 



V' 



(4.8) 



1-10 
0-110 


States of subsystems (i.e. individual photons) by the von Neumann's reduction (1.1), (4. 2), 
(4.3) are: 

1 
1 



(A 



aN) 



{Pa min) 



, or (pisn) = {hmin) = 7; 

1 2 



(4.9) 



V 
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By (2. 6), (4.4), (4. 5), for both subsystems it is equivalent to the replacement of the density 
matrix (4.8) by the reduced density matrix 

(Ps) (PSRn) = Pmin = ^1- (4.10) 

The expressions (4.9) and (4.10) show that each photon can be registered in any eigenstate 
with equal probability (1/2) but they do not determine the mutual correlation of results 
of simultaneous or consecutive measurements of both photons. The approximate nature 
of the von Neumann's reduction makes itself evident in the nonequality (4.8)7^(4.10). On 
the other hand, as photons in the pure state (4.8) are noninteracting, one might expect 
the exact expression of the system (4.8) through correlated subsystems in terms similar 
to (3.3), (3.4) to exist. With (4.6) and (4.7) for (4.8) the iterative procedure (3.3), (3.4) 
converges already in the first iteration for the arbitrary zeroth-order iteration and results 
in the following relations between the elements of correlated reduced density matrices of 
subsystems: ^22 = Ai; <^2i = "~/32i, in addition to usual relations of normalizing and of 
hermiticity: an = 1 — 0:22, cn2 = a2i- With these correlation conditions, let us write down 
the system of equations for calculating the elements 0:22, Oi2i from (2.1). This system is 
separated into three groups of equations: for diagonal elements, nondiagonal elements and 
equations of their relation: 

«22/522 ^ =: 022(1 - «22), «2l/5l2 ^ "1/2 =: -|tt21p, "22/^21 ^ =: Q;22(-«2l), 



Q;2i(1 - Q!22), 

Q;2ia;22, 

-(1 - Q;22)a2i- 



Q!22/5ll ^ 1/2 =: Q!22, Q!2l/?21 "S"^ =: - (q!2i)^, Q!21/322 ^ 

aii/322 ^ 1/2 =: (1 - q;22)^, Q;2i/3ii ^ =: 

aiiAi a22)oi22, Q;ii/32i ^ 

(4.11) 

The system (4.11) is obviously incompatible as an algebraic system of equation, but it is a 
compatible system of statistical equations that we have designated by "=:" instead of 
For example, the system of four equations from the first column of (4.11) is compatible, if 
it is granted that ^22 takes two values: 022 = 1 ^"^^ '^22 — with equal probability (1/2), 
and each equation is considered as a result of statistical averaging. In this way, for the first 
equation of the first column from (4.11) we have: = (1/2) q;22'*(1 ~ '^22'') + (^22^ ~ '^22'*) ■ 



For the second equation: 1/2 = (1/2) 



In solving the system (4.11) as a 



system of the statistical equations, we obtain within the arbitrary phase factor (exp(i^)): 

{Ps) = jlv^{[p%)®[p%)] 
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J_ id 

V2^ 



V2^ 



V 



J_„-i(6»+7r) 



+ 



J_„i(6l+7r) 



^e- 



J_ -i(e+7r) 



73' 



+ 



+ 



V %/2 



J_„-i(6(+7r/2) 



]_g-i(6»+37r/2) 



J_„i(6»+7r/2) 



J_„i(6»+37r/2) 



V2^ 



J_„-i(6(+37r/2) 
V ^/2^ 



J_„-j(e+7r/2) 
%/2^ 



J_pj(e+37r/2) 
V2^ 



J_ i(e+7r/2) \ 
V2^ 





(4.12) 



Consequently, the entangled pure EPR state is exactly represented as a sum of four hidden 
equiprobable (p^ = const =1/4) states of separated but correlated {p^aC^P^pc) subsystems. 
Each of summands presents by itself the distinct state of EPR pare. Each of p^*^, ff^Q cor- 
responds to the existence of the subsystem in one of eigenstates: the diagonal elements of 
density matrices are 1 and only. Note simultaneously the unusual nonzero values of non- 
diagonal elements. Nondiagonal elements of this type reflect not the usual superposition of 
eigenstates of one subsystem but the correlation in phase of the eigenstates of the correlated 
subsystems. 

The considered model of a pair of noninteracting two-level subsystems is the elementary 
example of the compound closed quantum system. The obvious next step of the generaliza- 
tion of the model is taking into account the interaction between subsystems. 



B. Dynamics of a pair of interacting two-level atoms 



The interaction of a pair of two-level subsystems was in detail considered in the theory 
of magnetic resonance, therefore it is convenient to take advantage of designations accepted 
there. The hamiltonian of a pair of identical spins- 1/2 in an external d-c magnetic fleld with 
a simple form of the spin-spin interaction is [14] : 



u+J d 
-J c 
c -J 
d -u + J 



(4.13) 
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where hcu is the Zeeman energy of a single spin- 1/2 in an external d-c magnetic field, pa- 
rameters J, c, d characterize the energy of the spin-spin interaction, which for definiteness is 
believed to be rather weak (a; > J, c, (i). Designations of basic states accepted in (4.1) are 
used, i.e. |2), |1) are eigenstates of a single spin-1/2 in the d-c magnetic field with higher 
and lower energies, respectively. 

The operator of evolution of quantum system (4.13) results in: 



exp ( T^^H 



( 



\ 



m 

U\t) 

C{nt) ^ f -S(Qi) 



^{id/n)e^'-^*S{VLt) 



T{id/^)e^'^*S{Vtt) 

^±iJt cos(ct) ^ie=^*"^* sin(ci) 



^ie^^*-^* sin(ci) e=^'"^* cos(ct) 















(4.14) 



where C{Qt) — cos(fii), S{Qt) — sm{Qt). Thus, the density matrix {p{t)), i.e. the dynamics 
of the system (4.13) at the given initial condition (p(0)) is known exactly. To describe the 
dynamics of each spin one should calculate correlated reduced density matrices of spin-1/2 
subsystems according to Sec. III. Let the initial state of the system (4.13) be: 







COS^ if 





— sin (f cos (f 

sin if cos ifi sin^ (p 





(4.15) 



For example, the initial condition of a radical pair born during photosynthesis is described 
in this form [15]. Similarly to the consideration of Subsec. IV A, this initial density matrix 
is presented as: 



p(0) ~ pRc(0) 



— < 



cos^ ip 



COS^ if 




+ sin if 




Ps, 
Pt, 



, I tanc^l > 1, 

, |tan(/7| < 1, , (4.16) 

tan(/7 = 1, 
tan(/7 = —1, 
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where 



(Pt) 



'^O 0^ 
110 
110 




1 

4 



1 -^e'^ 







71' 



75^^' 



+ 



+ 



+ 







J_pi(e+7r/2) 



J_ -i(0+7r/2) 



, J_p-i{e+T/2) 



J_pj(e+7r/2) \ 





J_„i(0+7r) 



J_„-i(6»+7r) 



1 



V2' 



1 -j(e+7r) 
V V2^ 



1 „j(6l+37r/2) 
V2^ 



1 „-i((9+37r/2) 
V V2^ 



1 „-i(6»+37r/2) 
V V2^ 



1 „i(6»+37r/2) 

v/2^ 





(4.17) 



Unlike (4.12) expression (4.16) is exact only for |tan(/?| = 1 (i.e. for the last two rows in 
(4.16)). 

As a result of correlated reduction (3. 3), (3. 4) the density matrix p{t) is represented at 
any time as follows: 



m) - 2 





Piv,t) icos(2<^)sin(2ct) -sin(293) 

-icos(2(^)sin(2ct) -sin(2(^) M{ip,t) 









{PRc{t) 



1 


1 




+ M{<p,t) 





C((^,i)<Q(,4.18) 

tan((^) — 1, 
tan((/?) = —1, 

where C{ip, t) = cos(2(^) cos(2ct), P{(p, t) = 1 + C{ip, t), M{ip, t) = 1 - C{(p, t). 

Thus, the dynamics of the closed system of two interacting spins-1/2 (4.13), according 
to the initial state (the value of (/?), may be of two types. For | tan(/?| = 1 the system is in 
the stable statistically defined state S<->(4.12) or T<->(4.17). For |tan</7| 7^ 1 there are the 
harmonic oscillations of probabilities to detect spins in basic states opposite in phase. At 
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tn — (2n + l)n/8c, {n = 0, 1, 2, . . .) for an arbitrary initial state these probabilities are the 
same and equal 1/2. 



V. DYNAMICS OF A TWO-LEVEL ATOM IN A RESONANT FIELD 



To demonstrate the decoherence phenomenon, the model problem of quantum optics on 
the dynamics of the inversion of a two-level atom interacting with a resonant electromagnetic 
field of a single- mode cavity (see, e.g., [2, 5]) is used traditionally. It is attractive due to 
the existence of the experimental setup of the cavity quantum electrodynamics entangling 
machine [5] and the existence of the exact solution of the theoretical Janes-Cummings model 
(JCM) [16, 17]. 

Let us trace basic distinctions in the dynamics of atom and field following from definitions 
(1.1)^(2.6) and (2.2), (2.3)^(3.3), (3.4). 

The hamiltonian of the resonant JCM looks hke [16, 17, 18, 19, 20, 21]: 

H = Ha + Hf + Haf = ^ {P22 - Ai) + ^ (a+a + aa+) + ^ (Aici - A2a+) , (5.1) 

where Ha, Hf, Haf are hamiltonians of the atom, field and atom-field interactions, cu is 
the resonant frequency of the field, Q is the constant of the atom-field interaction, Aj 
arc projection atomic operators representing square matrices of the second order, the ij- 
th element of which is equal to 1, others are equal to 0, i,j = 1,2 and the subscript 2 
corresponds to the level with the higher energy. 

It is convenient to write down the solution of the JCM (obtained in a general form in the 
Heisenberg picture in Ref. [19]) in the Schrodinger picture in the form of expression for the 
operator of evolution [20, 21]: 



> = exp ^—H = F22 exp ( ^lujtaa cos — ( aa 
C/t(i) \ n I \ J 9. y 



1/2 



+P11 exp (^^icuta'^ cos — (ji^a^ ^ 
^A2exp (^^iujta'^ cij exp(—i(p) sin ^ 

where exp(±i(^) = 



± Ai exp (^^iujtaa^^ exp{i(p) sin — ^a"*"a^ 

1/2 



1/2 



aa 



(5.2) 



> is the operator of a phase of a field [22]. Use of 



(a+a + l)"'/'a 
a+ (a+a + l)-'/' 

the operators exp(±i(^) and Aj (instead of more widespread Pauli operators) allows us to 
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write down the operator of evolution in the symmetric compact form (5.2), convenient for 
the further analysis. 

The elementary example (traditionally used in educational purposes [16, 17]) is the 
consideration of the dynamics of the inversion of the atom, if the initial condition was 
p(0) = Pa(0)p/(0) = P22|0)(0|, i.e. the atom was in the excited state |2) and the field was 
in the vacuum state. The known result, at once following from (5.2) and (1.1), consists 
in oscillations with the frequency Q, of the inversion of the atom or the probability of the 
existence of the photon in the cavity (vacuum Rubi oscillation [5]): 

PaNit) = P22 cos\flt/2) + Ai sm\flt/2), (5.3) 

or 

PfN{t) =sin2(Oi/2)|l)(l|+cos2(OV2)|0)(0|. (5.4) 

At t ^ 2-Kl/il, (/ = 0, 1, 2, . . .) the atom and the field are in a mixture state corresponding 
to decoherence of the observed subsystem. 

It is easy to notice that expressions (5.3) and (5.4) are incompatible for the system 
(2. 2), (2. 3). We search for the correlated reduced density operators of the atom and the field 
with (3. 3), (3. 4), using (5.3) and (5.4) as the zeroth-order iteration. As a result, we obtain: 

Pa{t) = C{t)p22 + S{t)Pn, Pf{t) = S{m{0\ + C{t)\l){l\, (5.5) 
where S{t) = 1 - C{t), s = sin(f]t/2), c = cos(l]t/2). 



C{t) = lim 



^2n 



.2n I o2n 



n— »oo Q'^^ -|- g 



1, [{Al - l)7r]/2Q < t < [(4/ + l)7r]/2Q, 
1/2, t = [{21 + l)TT]/2n, (5.6) 

0, [(4/ + l)7r]/21] < t < [{41 + 3)7r]/21^. 



According to (5.5) and (5.6), probabilities to find the atom and the field in their eigenstates 
are stepwise periodic functions of time with the period ir/fl. Practically all time the atom 
and the field are in pure states and only at the discrete moments of time they are in a 
mixture state equiprobable for two probable pure states. 

The expression (5.3) is employed at once to describe the experimentally observed vacuum 
Rubi oscillations [5, 23]. We call attention to the experimental "dead time" prior to the 
first-time photon emission in these experiments. The authors of [5, 23] included the effective 
interaction time taking into account the spatial variation of the atom-field coupling when 

19 



atom moves across the cavity. It seems us unreasonable for the initially vacuum field up 
to the first photon emission. All the space: the interior-space and the outer-space of the 
cavity, is uniform for the vacuum field. Thus all the atom way after the excitation must be 
included into consideration and not just the interior of the cavity. Then the experimental 
results [5, 23] contain the dead time when atom is in the excited state and there is not 
spontaneous emission that corresponds qualitatively to the time interval (0 7r/20) of 
(5.5),(5.6). 

The other experimental result corresponding qualitatively to (5. 5), (5. 6) consist in the 
observing of the omnidirectional superfiuorescense [26, 27] . There is the well-known exper- 
imental situation, when each atom of the system of many uncoupled two-level atoms for a 
short (in comparison with tt/Q) time interval is transferred into the excited state and the 
spontaneous emission of the inverted system of atoms is investigated without cavity. The 
pulse of radiation experimentally observed in this case after a time interval r is traditionally 
interpreted, after the authors of first observation [24], as a pulse of superradiation or super- 
fiuorescence [25]. And r is considered to be defined by the time of self-organization in the 
system of atoms at the expense of their interaction through a common radiation field [24]. 

According to (5.6), each atom emits a photon at the same time r = 7r/2f2 if at t = 
each atom was in the excited state, i.e. the pulsed emission of spontaneous radiation by 
the system of noninteracting atoms may be due to the correlation of each atom with the 
initially vacuum field. And for the gaseous spherical sample [26, 27] the directions of the 
spontaneous emission of photons are equiprobable. 

VI. CONCLUSION 

Results of the present work may be formulated as the following assertions: 

1. The identification of an operator of an observable of a subsystem of a closed quantum 
system with the appropriate extended operator of the independent observable subsys- 
tem {A' A) corresponds to the physical approximation, at which this obserable is 
calculated provided that the unobservable subsystem is in the steady state of minimum 
information (infinite temperature). 

2. The known algorithms of reduction correspond to the approximation of a given state 
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of the unobservable subsystem. 

3. The algorithm of correlated reduction that most corresponds to the ideal quantum 
measurement of observables of subsystems of the closed quantum sistem may be pre- 
sented as the self-congruent calculation of mutually correlated reduced density opera- 
tors of subsystems Pacit), Ppcit) by the successive approximation method (3. 3), (3. 4). 
The mean values of observebles of subsystems at the ideal quantum measurements are 
defined by expressions (^(^))^ = Sppac{t)A, (-B(t))^ = Spppc{t)B. 

4. The correlations in the dynamics of subsystems are demonstrated with examples of 
simple compound quantum systems: 1) the pair of the coupled two-level atoms (spins- 
1/2), and 2) the two- level atom in a single- mode resonant field. 

Acknowledgments 

The author acknoledge numerous discussions with V. N. Lisin. 



[1] D. N. Klyshko, Uspechi 168, 975 (1998). 
[2] M. B. Mensky, Uspechi 170, 631 (2000). 
[3] S. Ya. Kilin, Uspechi 169, 507 (1999). 

[4] I. Prigogine, The End of Certainty (The Free Press, NY-London- Toronto-Sidney-Singapore, 
1997). 

[5] J. M. Raimond, M. Brune and S. Haroche, Rev. Mod. Phys. 73, 565 (2001). 

[6] J. von Neumann, Mathematical Foundations of Quantum Mechanics (Princeton University 

Press, Princeton, 1955). 
[7] P. A. M. Dirac, The Principles of Quantum Mechanics (Clarendon Press, Oxford, 1958). 
[8] U. Fano, Rev. Mod. Phys. 29, 74 (1957). 

[9] K. Blum, Density Matrix. Theory and Applications (Plenum Press, New York and London, 
1981). 

[10] N. K. Solovarov, in Proceedings 2 of Intern. Conf. Geometrization of Physics V, edited by 
V. I. Bashkov, (Kazan State University, Kazan, Russia, 2001), p. 112. 



21 



[11] N. K. Solovarov, in Zavoisky Physical- Technical Institute, Annual 2001, edited by K. M. Sa- 

likhov, (Phystech Press, Kazan, 2002), p. 104. 
[12] N. Imoto, M. Ueda and T. Ogawa, Phys. Rev. A41, 4127 (1990). 
[13] T. Fukuo, T. Ogawa and K. Nakamura, Phys. Rev. A58, 3293 (1998). 

[14] A. Abragam and B. Bleaney, Electron Paramagnetic Resonance of Transition Ions (Clarendon 

Press, Oxford, 1970), v. 1, ch. 9.5. 
[15] Yu. E. Kandrashkin, K. M. Salikhov and D. Shtelik, Appl. Magn. Reson. 12, 141 (1997). 
[16] W. H. Louisell, Quantum statistical properties of radiation (John Wiley &; Sons, New York, 

1973). 

[17] L. Allen, J. H. Eberly, Optical Resonance and Two-level Atoms (John Wiley and Sons, NY- 
London- Sydney-Toronto, 1975). 
[18] H. I. Yoo and J. H. Eberly, Phys. Reports 118, 239 (1985). 
[19] J. C. Ackerhalt and K. Rz^zewski, Phys. Rev. A12, 2549 (1975). 
[20] N. K. Solovarov, Optica Acta 27, 393 (1980). 
[21] D. A. Demidov and N. K. Solovarov, Laser Physics 5, 997 (1995). 
[22] R. Loudon, The Quantum Theory of Light (Clarendon Press, Oxford, 1973). 
[23] M. Brune, E. Hagley, J. Dreyer et al, Phys. Rev. Lett. 77, 4887 (1996). 
[24] N. Skribanovitz, I. P. Herman, J. C. MacGillivray et al., Phys. Rev. Lett. 30, 309 (1973). 
[25] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[26] A. L Lvovsky, S. R. Hartmann and F. Moshary, Phys. Rev. Lett. 82, 4420 (1999). 
[27] A. L Lvovsky, S. R. Hartmann and F. Moshary, Phys. Rev. Lett. 89, 263602 (2002). 



22 



